A new synchronization mechanism via Turing-hke microscopic structures 

for CO oxidation on Pt(llO) 
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We discuss an alternative to the traditional gas-phase coupling approach in order to explain 
synchronized global oscillations in CO oxidation on Pt(llO). We use a minimalist microscopic 
model which includes structural Ft surface reconstruction via front propagation, and large diffusion 
rates for CO. The synchronization mechanism is associated with the formation of a Turing-like 
structure of the substrate. By using large parallel microscopic simulations we derive a scaling laws 
which allow us to extrapolate to realistic diffusion rates, pattern size, and oscillation periods. 

PACS numbers: 82.65.+r; 82.20.Wt; 02.70.Tt; 82.40.Np; 89.75.Da 
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Oscillations appear in experiments of CO oxidation on 
Pt(llO) in a very narrow parameter interval Addi- 
tionally, it is also possible to observe pattern formation; 
about T ^ 500K a crossover exist from spirals and fronts 
at low T to standing waves and chemical turbulence at 
high T . Experiments using fast scanning tunneling mi- 
croscopy and field ion microscopy have clearly shown that 
fast kinetic processes are typically accompanied by the 
appearance of microscopic structures ^, ^ . 

Although the diffusion of CO is fast, it does not provide 
a synchronization mechanism for the global oscillations, 
it is thought that the mechanism for global synchroniza- 
tion is coupling via the gas-phase ||^. This mechanism 
has been used widely in classical reaction-diffusion (RD) 
models. However, experiments have shown that particle 
dynamics on surfaces can be highly correlated j^] , some- 
thing which is not taken into account in RD models. This 
leads to differences between RD and experiments which 
can only be removed by including phenomenological cor- 
rections. An alternative mechanism for global synchro- 
nization of oscillation has been suggested based on spon- 
taneous phase nucleation Q . This nucleation results in 
a random creation of dynamic defects and leads to global 
synchronization via stochastic resonance. 

We present in this paper Monte Carlo (MC) simula- 
tions of a microscopic model for CO oxidation on Pt(llO). 
We show that our model can produce global oscillations 
via a new mechanism showing microscopic structures 
without gas-phase coupling, coupling via stochastic reso- 
nance, or lateral interactions. The microscopic structures 
appear in the substrate, but they also support pattern 
formation in the adlayer on a mesoscopic length scale. 
In addition we provide a connection between these MC 
simulations and experimental results, via scaling laws. 

MC simulations provide a systematic approach to in- 
clude microscopic effects MC simulates directly the 
chemical reaction steps in a model, without averaging out 
particle correlations. Simulations of our model were per- 
formed by using a Cellular Automaton (CA) which has 



been shown to be equivalent to MC A limitation of 
MC simulations has been the inability to deal with exper- 
imental system sizes and realistic diffusion coefficients. 
However, by using large parallel simulations as described 
in pO| we can now estimate results for experimental con- 
ditions. Using this parallel code we are able to simu- 
late system sizes of about I /im with diffusion coefficients 
of about 10~^''cm^sec~^. This approaches experimental 
conditions closer than any other previous MC simulation. 
Experimental values are for system sizes 10^ jiva or more 
and diffusion coefficient 10~^cm^sec~^. There is a gap 
between experimental and simulated rate constants, but, 
as we can vary the system size and the diffusion rate over 
a large range, we can derive scaling laws, which we can 
use to extrapolate the behavior predicted by our model 
to realistic system sizes and diffusion rates. 



Our model has already been used in several studies 
1^, |ll|, O2 adsorbs dissociatively onto two nearest 

neighbor sites with rate constant (1 — y)sy^ with x = Oi^Pi 
where a denotes the 1x2 phase and (3 denotes the unre- 
constructed 1x1 phase of Pt(llO) ||l|. The experimental 
value for the ratio of the sticking coefficients of O2 on the 
two phases is Sa '■ sp ~ 0.5 : 1. CO is able to adsorb onto 
a free surface site with rate constant y and it desorbs 
from the surface with rate constant k. Both reactions 
are independent of the surface phase to which the site be- 
longs. In addition CO is able to diffuse via hopping onto 
a vacant nearest neighbor site with rate constant D. The 
C0+0-^C02 reaction occurs with rate constant R, when 
CO and O are nearest neighbors sites. CO2 desorbs im- 
mediately forming two vacant sites. The a ^ [3 surface 
phase transition is modeled as a front propagation with 
rate constant V . For two nearest neighbor surface sites in 
the state a(3 the transition a(3 aa {a/S — > /?/?) occurs 
if none (at least one) of these two sites is occupied by CO. 
Summarizing the above transition definitions written in 



the more usual form of reaction equations gives: 



CO{g) + S^ ^ CO(a), 
02(g) + 25"^ 20(a), 
02(g) + 25/5 - 
CO(a) + S^~- 
CO(a) + 0(a) 
3°" ^ S*^, 



20(a), 
S''^ + CO(a), 
-^C02(g) + 2^x, 



where S stands for a vacant adsorption site, x stands for 
either a or /3, and (a) and (g) for a particle adsorbed on 
the surface or in the gas phase, respectively. For addi- 
tional details see Ref. . Our simulations a use special 
normalization of the rates: the adsorption rate for CO is 
related to the partial pressures by 
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Pco + Po^ 



(1) 



Comparing our rate constants with typical experimental 
values used in RD models jl^, |l^ we can estimate the 
unit of time in our simulation to be about ~ 10~^s. 
All our rate constants are expressed in that unit of time. 
Lengths are expressed in units of the unit cell parameter 
a. The diffusion rate D is related to the experimental 
diffusion coefficient Da by Da = a^D/zto, where z = 4 
is the coordination number. 

The quality of the oscillations as a function of y has 
a typical resonance form. Well-defined oscillations are 
found near y — 0.494. Under certain condition, which 
we discuss below, these oscillations are global. Moving 
away from y = 0.494 we see first that the synchronization 
decreases and then a disappearance of the oscillations. 
We restrict ourselves to y = 0.494 (for a fixed values 
k = 0.1 and = 1) in our simulations. 

We have found that two situations are possible depend- 
ing the system size, the diffusion rate, and on the initial 
conditions (IC). If the system is small or the diffusion 
is fast we find global oscillations. If the diffusion is too 
slow for a given system size the IC become important. 
We can again have global oscillations or the formation of 
patterns in the form of moving fronts and spirals. 

To illustrate the synchronization mechanism Fig. |l| 
shows global oscillations. When the system is covered 
almost fully by CO the (3 islands in the a background 
grow (snapshots 1 and 2). When the size of the (3 islands 
is large enough the rate of adsorption of O2 becomes 
larger than the adsorption of CO on the /3 zones (snap- 
shots 3 and 4). At the same time the CO in the a zones 
is converted to CO2 near to the borders with the (3 zones 
(snapshot 4) . When the total coverage of oxygen becomes 
larger than CO, then the direction of the island grow is 
reversed (snapshots 4, 5, and 6) and a islands grow in the 
P background. When the size of the /? zones is not large 
enough to keep the oxygen coverage larger than the CO 
coverage (snapshot 7) , then both phases are become cov- 
ered by CO, and the cycle start again (snapshot 8). The 
key point in this cycle is the existence of a critical size 
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FIG. 1: Oscillations for D = 4000 and L = 1024 on the cov- 
erage CO (solid), oxygen (dashed), and a (dotted). Sequence 
of temporal snapshots correspond to the points in the upper 
figure. Each snapshot has two parts: In the left part we plot 
the chemical species; CO particles are grey and O particles 
are white, and empty sites are black. The right part shows 
the structure of the surface; a phase sites are black, and /3 
phase sites are white. 



of the P islands where oxygen becomes more stable than 
CO and a critical size of a islands where CO is more sta- 
ble than oxygen. This corresponds to a phase transition 
produced by a varying size of these islands. This phase 
transition is driven by diffusion and the critical size of 
the islands depends on the diffusion rate as a/D. 

We found that the most important feature of the struc- 
tures of the substrate is the typical distance between the 
islands; i.e., the correlation length Ic of these structures. 
For Ic we use the distance to the first maximum in the ra- 
dially averaged correlation function of the substrate. The 
shape of the islands self-organizes through successive os- 
cillations so that they cover the whole system with same 
sized islands separated by similar distances and small 
differences in size are removed. We note that our sim- 
ple model of surface reconstruction by border propaga- 
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FIG. 2: Global oscillations and pattern formation with 
D — 250. Sections of the upper-left corner with L = 8192, 
4096, 1024, and 256 are shown on the left side. The sections 
correspond to the dot in the temporal plot at the top. It 
shows the same information as Fig. 0, but it also includes the 
rate of CO2 formation R. On the right side we have a wave 
front, a target, a spiral, and turbulence (L = 2048), which 
can be obtained with different initial conditions. 

tion dynamic is sufficient to produce this effect. It yields 
synchronization at large scales even with relatively small 
diffusion rates. In Fig. ^ the diffusion length f = V DT, 
with T the period of the oscillations, is larger than the 
system size L. The synchronization is therefore trivial. 
In Fig. H, however, a simulation is shown where the dif- 
fusion length is much smaller than the lattice size, but 
larger than the correlation length. So the diffusion syn- 
chronizes the oscillations on neighboring islands, and the 
self-organization leads to global oscillations. 

In cases where the synchronization is even better, as 
with large D values, islands tend to be grouped form- 
ing rolls which are parallel (see Fig. |l|). The structures 
formed in this way resemble labyrinthic Turing struc- 
tures. Moreover a percolating connection between islands 
can be observed. The full analytical analysis of this be- 
havior is far from trivial, however the difference of the 
mobility for CO and O, could be sufficient for get a Tur- 



ing instability in this model. 

When the diffusion is slow pattern formation can also 
occur as on the right of Fig. ||. The distance between the 
centers of a and f3 islands is the same as for the global 
oscillations, but the size of the islands differ for differ- 
ent parts of the system. These differences will not be 
removed as in the fully synchronized cases. As a conse- 
quence a delay or detuning will appear in the oscillations 
for different parts of the surface, which can create fronts, 
target pattern patterns, spirals, and turbulence. A cross- 
section of one of these fronts shows a delay between the 
pulse of the adsorbates and the pulse of the substrate 
phases, which determines the direction of propagation. 
This is a common feature of excitable, bistable systems, 
such as the activator-inhibitor models. The most remark- 
able feature is that these fronts appear on top of fixed 
Turing-like structures which are present even in these in- 
completely synchronized cases. So we have within our 
model two different length scales for pattern formation; 
excitable dynamics for fronts and spirals on a mesoscopic 
scale (Hopf instability, 10), and quasi- fixed Turing- like 
microscopic structure. 

A scaling analysis for the Turing structures and the 
oscillations for the fully synchronized case is presented 
in Fig. H We show the scaling behavior for the pe- 
riod of oscillations T and the correlation length Ic as a 
function of diffusion rate D. Numerically we found that 
the relations T{D) - D" and lc{D) - holds with 
ly = 0.4849 ± 0.016 and = 0.5069 ± 0.01. For large D 
values a large system size is required to compute lc{D) 
precisely. For instance for D = 8000 we use L = 8192. 
Additionally large diffusion rates means that most of the 
time is spend moving particles; on 8000 diffusion steps 
only one chemical reaction happens. Also the simula- 
tion time increase because the oscillation period increases 
with D. 

A relationship T{D) ^ lc{D) ^ D^/^ seems quite rea- 
sonable, as our model has fast diffusion. In fact the scal- 
ing law for Ic is the consequence of scaling law for T . 
The width of the Turing structure oscillates on a length 
of order Ic due to the phase propagation with rate V . 
The period is related to these by Ic ~ VT . If the velocity 

V is constant, then Ic ~ T. Fig. ^ was computed for a 
fixed value of V . Using a larger value for the velocity pa- 
rameter we can obtain larger correlation lengths Ic and 
smaller periods T, but the scaling remains the same. We 
have found for the Turing-synchronized oscillations that 
the relation T ^ V^^/"^ also applies as in |^, |l^. So 
we have an important scaling law in form T — c^J D/V, 
where c ~ 10. An estimate for the oscillation period 
for realistic diffusion rates Da ~ 10~^ cm^sec"^ we have 
D 10^. The experimental oscillations period is r = igT 
with to 10~^ sec and r = 10 sec. The value V predicted 
by our model for that value of t is then of the order of 

V - 102. 

We propose that the minimal condition to have global 
oscillations is that the diffusion length ^ is at least of the 
order of the correlation length Ic'- S, = Ic {S, = V DT and 
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FIG. 3: Scaling behavior for the period of oscillations T{D) 
and the correlation length lc{D) (units). 

Ic = VT, c ^ 10). Critical values for the diffusion rate 
and the correlation length are then 

D = cV^ (2) 

For V — 1 we get the order of the parameters used in 
Fig. I (left side), D ~ 10^ and Ic ^ 10^ where global 
oscillations are possible but infrequently. 

Diffusion is a thermally activated process so we can 
compare our proposal, Eqs. (H), with the experimental 
crossover at T ~ 500; for lower temperatures where fronts 
and spirals appear, we have lower diffusion rates, f < 



Ic and the synchronization mechanism is not stable and 
produces patterns similar to on the right of Fig. ||. On the 
other hand for large temperatures where standing waves 
and chemical turbulence appear, we have large diffusion 
rates, f > Zc, and we get better synchronization, as in 
Fig. hi. For ^ = /c we get the crossover, Eqs. d) holds 
and we get the minimum criteria for synchronization, as 
in Fig. g on the left. 

Using V — 10^ in Eqs. (H) we get for the critical value 
of diffusion rate D ~ 10* , which corresponds to the ex- 
perimental value of Da 10~^cm^sec~^. This confirms 
our crossover idea. For correlation length 1^ ^ 10^ (in 
units of cell parameter a) we have a size ~ 10^ /im which 
is of the order of magnitude of the standing waves ob- 
served in experiments |]l^. Near crossover it is possible 
to interpret V also as the velocity v of the fronts and 
spirals. For V — 10^ we have v = aV/to ~ 10~^cm s~^, 
in agreements with experiments [T^ . 

To summarize we present in this paper a new mecha- 
nism for synchronization of global oscillations based on 
microscopic Turing-like structures in the reconstruction 
of the surface. In cases with incomplete synchronized 
oscillations the adlayer has a mesoscopic second charac- 
teristic length. This is the first actual demonstration of a 
double length scale. Scaling laws are analyzed to connect 
the model to experimental parameter values. 

This work was supported by the Nederlandse Organ- 
isatie voor Wetenschapperlijk Onderzoek (NWO), and 
the EC Excellence Center of Advanced Material Research 
and Technology (contract N lCAl-CT-2080-7007). We 
would like to thank the National Research School Com- 
bination Catalysis (NRSCC) for computational facilities. 



[1] R. Imbihl, G. Ertl, Chem. Rev. 95, 697 (1995). 

[2] J. Wintterlin, Chaos, 12, 108 (2002). 

[3] S. Volkening, K. Bediirftig, K. Jacobi, J. Wintterlin, 

G. Ertl, Phys. Rev. Lett. 83, 2672 (1999). 
[4] J. Wintterlin, S. Volkening, T.V.W. Janssens, T. Zam- 

belli, G. Ertl, Science 278, 1931 (1997). 
[5] H. Levine, X. Zou, Phys. Rev. Lett. 69, 204 (1992). 
[6] O. Kortliike, V.N. Kuzovkov, W. von Niessen, Phys. Rev. 

Lett. 83, 3089 (1999). 
[7] O. Kortliike, V.N. Kuzovkov, W. von Niessen, Phys. Rev. 

E66, 036139 (2002). 
[8] J.J. Lukkien, J.RL Segers, P.A.J. Hilbers, R.J. Gelten, 

A.RJ. Jansen, Phys. Rev. E 58, 2598 (1998). 
[9] O. Kortliike, J. Phys. A 31, 9185 (1998). 
[10] R. Salazar, A.P.J. Jansen, V.N. Kuzovkov, 



nlin.CG/020705c . 
[11] O. Kortliike, V.N. Kuzovkov, W. von Niessen, Phys. Rev. 

Lett. 81, 2164 (1998). 
[12] V.N. Kuzovkov, O. Kortliike, W. von Niessen, Phys. Rev. 

Lett. 83, 1636 (1999). 
[13] O. Kortliike, V.N. Kuzovkov, W. von Niessen, J. Chem. 

Phys 110, 11523 (1999). 
[14] M. Bar, M. Hildebrand, M. Eiswirth, M. Falcke, H. Engel, 

M. Neufeld, Chaos 4, 499 (1994). 
[15] A. von Oertzen, H.H. Rotermund, A.S. Mikhailov, 

G. Ertl, J. Phys. Chem. B 104, 3155 (2000). 
[16] M. Falke, M. Bar, H. Engel, M. Eiswirth, J. Chem. Phys 

97, 4555 (1992). 



